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Pulsars are rapidly-rotating, highly-magnetized neutron stars emitting radiation across the elec- 
tromagnetic spectrum. Although there are more than 1800 known radio pulsars, until recently, 
only seven were observed to pulse in gamma rays and these were all discovered at other wave- 
lengths. The Fermi Large Area Telescope makes it possible to pinpoint neutron stars through 
their gamma-ray pulsations. We report the detection of 16 gamma-ray pulsars in blind fre- 
quency searches using the LAT. Most of these pulsars are coincident with previously unidenti- 
fied gamma-ray sources, and many are associated with supernova remnants. Direct detection of 
gamma-ray pulsars enables studies of emission mechanisms, population statistics and the ener- 
getics of pulsar wind nebulae and supernova remnants. 

A wide variety of astrophysical phenomena, such as black holes, active galactic nuclei, gamma-ray bursts, and 
pulsars are known to produce photons exceeding many mega electron volts. Detection and accurate localization of 
sources at these energies is challenging both because of the low fluxes involved and the limitations of the detection 
techniques. The sky above 100 MeV was surveyed more than 30 years ago by the COS-B satellite (7), and more 
recently by EGRET (2), on the Compton Gamma-Ray Observatory. One of the main legacies of EGRET was 
the detection of ^300 gamma-ray sources, many of which have remained unidentified, despite searches at a wide 
variety of wavelengths (3). Many EGRET (and COS-B) unidentified sources are thought to be of Galactic origin, 
because of their lack of variability and concentration along the Galactic plane. A large fraction of these have been 
suspected to be pulsars [e.g. (4-6)] despite deep radio and X-ray searches often failing to uncover pulsed emission, 
even when the gamma-ray sources were coincident with supernova remnants (SNR) or pulsar wind nebulae (PWN). 
The lack of radio pulsations has usually been explained as the narrow radio beams missing the line of sight toward 
the Earth (7). We refer to such pulsars as "radio-quiet"; even though they may emit radio waves, these cannot be 
detected at Earth. Before the launch of Fermi, Geminga (8) was the only known radio-quiet gamma-ray pulsar. 
Current models of pulsar gamma-ray emission predict that gamma-ray beams are much wider than radio beams 
(9), thus suggesting that there may be a large population of radio-quiet gamma-ray pulsars. 

Soon after launch on 1 1 June 2008, the Large Area Telescope (LAT) on Fermi began surveying the sky at 
energies above 20 MeV. A companion paper describes the LAT detection of a population of gamma-ray millisecond 
pulsars (MSPs) (10). Here we report the detection of 16 pulsars found in blind frequency searches using the LAT. 
Previously, gamma-ray pulsars had been detected only by using a radio (or in the case of Geminga, an X-ray) 
ephemeris. 

Observations and Data Analysis. The LAT is a high-energy gamma-ray telescope sensitive to photon energies 
from 20 MeV to over 300 GeV, featuring a solid-state silicon tracker, a cesium-iodide calorimeter and an anti- 
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coincidence detector (77). Gamma-ray events recorded in the LAT have time stamps derived from a GPS clock on 
the Fermi satellite with an accuracy of < 1/is (72). Fermi operates in continuous sky survey mode, covering the 
entire sky every 3 hours. Compared to EGRET, Fermi has a larger effective area (9,500 cm 2 at normal incidence), 
larger field of view (2.4 sr), more efficient use of time on orbit for photon collection, and a finer point spread 
function (5° at 100 MeV, 0.8° at 1 GeV). The first three factors result in more rapid photon accumulation and the 
fourth increases the signal to noise ratio by improving the background rejection. 

Even with the improvements of the LAT, gamma-ray data remain extremely sparse. The brightest steady 
gamma-ray source in the sky, the Vela pulsar, spins 1 1 times per second; however, it provides fewer than one 
hundred (> 30 MeV) photons during every two orbits of the LAT (73). As a result, detection of gamma-ray 
pulsations from sources with more typical brightness requires weeks or months of data. A standard technique for 
finding a periodic signal in a data set is a fast Fourier transform (FFT). A fully coherent FFT becomes memory 
intensive because the number of frequency bins in the FFT increases with the length of the observational time, as 
A^bins = 2T/ max , where T is the duration of the observation and / max is the maximum frequency. Furthermore, 
pulsars gradually spin down as they radiate away energy, requiring the computation of many tens of thousands 
of FFTs to scan a realistic frequency and frequency derivative (/ and /) parameter space, making FFT searches 
computationally intensive, if not prohibitive [e.g. (74)]. We thus used the time-differencing technique (75) and a 
fixed window of T = 2 19 s (<~1 week). This dramatically reduces the number of FFT bins and frequency derivative 
trials. This method was successfully applied to EGRET data (76). 

We applied the time-difference technique to photons from a small region of interest around selected target 
positions. We used two sets of target positions. The first was a list of ^100 promising locations derived from 
multiwavelength observations. This list was selected based on location in the Galactic plane, spectra, lack of long- 
term variability, and presence of an SNR, PWN, or a likely neutron star identified from X-ray observations. Many 
of these locations were unidentified gamma-ray sources from EGRET. The second list, compiled after the LAT 
started collecting data, included locations of gamma-ray sources from a preliminary version of the LAT catalog. 
After eliminating sources with probable active galactic nuclei associations, we ended up with <~ 200 locations, of 
which 36 were consistent with previously unidentified EGRET sources. 

We analyzed data collected from sky survey observations beginning on 4 August 2008 and ending on 25 
December 2008. We applied a bary center correction using the JPL DE405 Solar System ephemeris (77), assuming 
that the source position was either a likely candidate counterpart, found through multiwavelengths studies, or our 
best estimate based on LAT data. This translation into the nearly inertial reference frame of the Solar System 
Barycenter (SSB) corrects for classical and relativistic effects of the motion of the Earth. We accepted all highest 
quality ["Diffuse" class (77)] photons falling within a 0.8° radius of our source location with energies greater than 
300 MeV. For each target, we searched a set of trial frequency derivatives from zero to the spin-down of the young 
Crab pulsar (/ = —3.7 x 10~ 10 Hz s _1 ). For each trial, we corrected the photon arrival times using the assumed 
frequency derivative, calculated the time differences and computed the FFT sampled with a Nyquist frequency of 
64 Hz, and searched the spectrum for significant peaks. 

We refined candidate signals by performing an epoch-folding search over a narrow region of frequency and 
frequency-derivative space using PRESTO (18). Using the parameters from the time differencing search as an 
initial timing model, we split the data into intervals (^2 — 4 weeks, depending on signal to noise), and measured 
pulse times of arrival (TOAs) in each segment. These TOAs were fit to a timing model using the Tempo2 pulsar 
timing package (79). We selected the reference epoch for all timing models to be MJD 54754.0, roughly in the 
middle of the observation interval. This choice minimized the covariances between the timing model parameters. 
We identified 16 pulsars within the first 4 months of data, and confirmed these using an independent data set of at 
least 1 additional month (Figs. 1 and 2; Tables 1 and 2). 

Six of the 16 pulsars were found using the positions of well-localized counterparts at other wavelengths, 
whereas 10 were found using LAT source positions. The LAT positions were uncertain by several arcminutes, 
and in many cases the initial timing models showed very large non-white residuals indicative of a significant dis- 
tance between the true position and the position being used to barycenter the data. We refined the positions in 
several ways (see below). The number of digits used in the declination of the names assigned to the pulsars reflects 
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our confidence in their localization. The discovery frequency may not be the true pulsar rotation frequency, but 
rather a harmonic or a subharmonic, especially for gamma-ray pulse profiles with two peaks separated by ~180 
degrees in phase. We tested for this effect, but for certain cases (e.g. J0357+32) the results remain inconclusive, 
because of the low photon counts. 

The gamma-ray pulsars. The 16 pulsars are representative of the highest spin-down luminosity portion of the 
general population, which includes 1800 pulsars in the ATNF database (20), including the 6 gamma-ray pulsars 
detected with EGRET, as well as 8 gamma-ray MSPs detected by Fermi-LAT (10). 

Thirteen LAT pulsars are associated with unidentified EGRET sources. Indeed, 15 out of 36 EGRET uniden- 
tified sources searched showed pulsations, but two of them (PSR J1028-5819 and PSR J2021+3651) are known 
radio pulsars (21, 22). Three pulsars (J0357+32, J1459— 60, and J2238+59) correspond to newly seen LAT sources. 
This result suggests that the blind searches are flux limited; although many LAT unidentified sources also could 
be pulsars, most have too low a flux for a pulsation to have been detected in our available data. Five pulsars are 
likely associated with PWN and/or SNRs and an additional one, J 1836+5925, is associated with a known iso- 
lated neutron star. J0007+7303 (23), which was long suspected of being a pulsar because of its clear association 
with SNR CTA 1 containing a PWN. J1418-6058 is in the complex Kookaburra region of the Galactic plane and 
is likely associated with PWN G313. 3+0.1, the "Rabbit" (24). J1809-2332, likely powering the "Taz" PWN, is 
also possibly associated with a recently-discovered mixed-morphology type SNR, G7. 5-1. 7(25) . J1826-1256 is 
probably powering the "Eel" PWN (26), whereas J202 1+4026 reveals that a pulsar is present, as long suspected, 
inside the gamma Cygni SNR (27). Two more pulsars present plausible associations with SNR: J0633+0632, in the 
fairly old (~ 30 kyr) complex Monoceros Loop SNR (G205.5+0.5) at a distance of - 1.5 kpc (28). J1907+0601 
is ~ 0.3° from SNR G40.5-0.5, but additional evidence is needed before an association can be demonstrated. 
J 1836+5925 is the long sought pulsar powering 3EG 1835+5918, the brightest and most accurately positioned of 
the unidentified EGRET sources (29). Our pulsar coincides with RX J 1836+5925, an isolated neutron star for 
which X-ray pulsations were extensively searched, but not found (30). Deep radio pulsation searches with the 
NRAO Green Bank Telescope (GBT) placed an upper limit at 1.4 GHz on flux density of 7 /iJy for P>10 ms (30). 
The X-ray counterpart of J 1835+59 18 was already suspected to be similar to Geminga, just a bit older and slightly 
farther away (within ^800 pc (31)). Using ~ 130 ks of publicly available Chandra HRC X-ray observations of RX 
J 1836.2+5925 (the same data searched by (30)), we searched for X-ray pulsations using our gamma-ray ephemeris, 
but found none. In addition to J1836+5925, three other pulsars are located off the plane, at |6| > 3° (J0007+7303, 
J0357+32, J1741-2054), whereas the remaining 12 are located in the Galactic plane. 

We imaged the error boxes of nine pulsars with specifically targeted short (<~5 ks) Swift observations, and 
in four cases identified X-ray point sources within the LAT error circle that may be the X-ray counterparts 
(J0633+0632, J1741-2054, J1813-1246, and J1958+2846). In one other case, J2032+4127, Chandra observations 
(32) have revealed a number of X-ray point sources in the region. The best timing solution was found assuming the 
position of MT91 221, but this cannot be the true X-ray counterpart as it is known to be a B star. The region around 
J1418-6058 was mapped in detail by Chandra (24), and a promising point source, denoted Rl, is consistent with 
the LAT source position and likely the X-ray counterpart to the pulsar. Chandra observations of 3EG 2020+4017 
show several possible counterparts to J202 1+4026 of which S21 (33) has the best timing solution and is also most 
consistent with the LAT source position. Other pulsars were better localized using grid searches (J 1459-60 and 
J1732-31) or on-pulse photon analyses (J0357+32, J1907+0601, and J2238+59). In our grid searches, we scanned 
25 equally-spaced locations overlapping the 95% LAT error circle, and looked for the position with the best timing 
solution. In on-pulse photon analyses, we used only the on-pulse photons to improve the source position. 

Five pulsars suffer from large (~ 0.1°) positional uncertainties, which degrade the frequency derivative in the 
timing solution. The problem is particularly acute in J0357+32 (Table 2), and the ~ 5' uncertainty in its position 
results in an uncertainty in the true value of the frequency derivative of ~ 6 x 10~ 14 . Any parameters derived from 
its timing solution are not reliable. 
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Implications. Data from EGRET, COS-B, and SAS-2 demonstrated that rotation-powered pulsars are persistent 
sources of pulsed gamma rays. EGRET results also showed that the luminosity of detected pulsars over all wave- 
lengths was dominated by the gamma-ray part of the spectrum. EGRET further established a list of persistent 
sources at low and intermediate Galactic latitudes that could not be demonstrated to show pulsations. The Fermi- 
LAT data confirms most of these sources, and shows that many of these persistent Galactic gamma-ray sources are 
pulsars. Of 36 LAT sources with EGRET counterparts that were searched, 15 show pulsations (13 from previously 
unknown pulsars). 

Until the launch of Fermi, the only known radio-quiet pulsar was Geminga, at a distance of ^250 pc (34), and 
with a spin-down luminosity of 3.2 xlO 34 ergs -1 and a gamma-ray efficiency of a few percent (for a 1 srbeam)or 
50% (for isotropic emission) (8). It was already apparent ~ 30 years ago, based on COS-B observations (1) that the 
average unidentified gamma-ray source should have a luminosity greater than 10 35 erg s _1 to be consistent with 
the very narrow distribution of low latitude gamma-ray sources (4, 35). Later studies with EGRET confirmed this 
(5), showing that most unidentified sources were at ~1.2-6 kpc, and had luminosities of 0.7 x 10 35 to 16.7 x 10 35 
erg s _1 . Thus, most unidentified sources could not be like Geminga [e.g. (8)]. The detection, by the LAT, of mostly 
high-luminosity, low Galactic latitude, and probably relatively distant pulsars confirms the predictions that the 
unidentified gamma-ray sources could not all be nearby, low-luminosity pulsars such as Geminga. 

This new population of gamma-ray selected pulsars helps reveal the geometry of emission from rotation- 
powered pulsars. A broad gamma-ray beam is required for at least part of our population of pulsars, as we see the 
high-energy beam but not the relatively narrow radio beam. Thus, this population favors pulsar emission models 
in which the high-energy radiation occurs in the outer magnetosphere, nearer to the light cylinder. Outer gap 
models [e.g. (6)], where particles are accelerated and radiate in vacuum gaps in the outer magnetosphere, and slot 
gap models [e.g. (36)], where particles are accelerated and radiate at all altitudes along the open magnetic field 
boundary, both can reproduce the wide, double-peaked light curves we observed (9). 
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Table 1: Names and locations of the gamma-ray pulsars discovered in our blind searches. The "Assumed Counter- 
part" column is the X-ray source that provided the position for the timing model, which in some cases may not be 



the true counterpart. 



LAT PSR 


LAT Source OFGL 


EGRET Source 


Assumed Counterpart 


RA(°) 


Dec (°) 


m 


b(°) 


J0007+7303 


J0007.4+7303 


3EG J00 10+7309 


RX J0007.0+7303 in Gl 19.5+10.3 (37) 


1.7565 


73.0523 


119.7 


10.5 


J0357+32 


J0357.5+3205 






59.445 


32.105 


162.7 


-16.0 


J0633+0632 


J0633.5+0634 


3EG J063 1+0642 


Swift J063344+0632.4 


98.4333 


6.5403 


205.0 


-1.0 


J1418-6058 


J1418.8-6058 


3EGJ1420-6038 


Rl inG313.3+0.1 (24) 


214.6779 


-60.9675 


313.3 


0.1 


J1459-60 


J 1459.4-6056 






224.874 


-60.878 


317.9 


-1.8 


J1732-31 


J1732.8-3135 


3EGJ1734-3232 




263.169 


-31.610 


356.2 


0.9 


J1741-2054 


J1742.1-2054 


3EGJ1741-2050 


Swift J174157-2054.1 


265.4908 


-20.9033 


6.4 


4.9 


J1809-2332 


J1809.5-2331 


3EG J 1809-2328 


CXOU J180950.2-233223 in G7.4-2.0 (38) 


272.4592 


-23.5397 


7.4 


-2.0 


J1813-1246 


J1813.5-1248 


GeV J1814-1228 


Swift J181323-1246.0 


273.3475 


-12.7668 


17.2 


2.4 


J1826-1256 


J1825.9-1256 


3EGJ1826-1302 


AXJ1 826. 1-1257 in G18.5-0.4 (26) 


276.5356 


-12.9429 


18.5 


-0.4 


J 1836+5925 


J 1836.2+5924 


3EG J 1835+59 18 


RX J1836. 2+5925 (30) 


279.0570 


59.4250 


88.9 


25.0 


J 1907+06 


J1907.5+0602 


GeV J 1907+0557 




286.965 


6.022 


40.2 


-0.9 


J1958+2846 


J1958.1+2848 


3EGJ1958+2909 


Swift J195846.1+2846.0 


299.6921 


28.7674 


65.9 


-0.2 


J2021+4026 


J202 1.5+4026 


3EG J2020+4017 


S21 inG78.2+2.1 (33) 


305.3777 


40.4462 


78.2 


2.1 


J2032+4127 


J2032.2+4122 


3EG J2033+4118 


MT91 221 in Cygnus OB2 (32) 


308.0612 


41.4610 


80.2 


1.0 


J2238+59 








339.561 


59.080 


106.5 


0.5 



Table 2: Rotational ephemerides for the new pulsars. For all timing solutions, the reference epoch is MJD 54754. 
The given frequency and frequency derivative are barycentric. The numbers in parentheses indicate the error in 
the last decimal digit(s). For LAT PSR J0357+32, marked with (|), the measured frequency derivative, and thus 
the derived parameters, may be significantly affected by positional error (see text).We also list the number of 
photons, n 7 , obtained with the cuts used in this work (see text) over the 5 month observational period. The (1-100 
GeV) photon number flux, F35, is given, as presented in the Fermi Bright Source List (39), except for LAT PSR 
J2238+59, which was not in the Bright Source List but whose flux was derived from the same data set and analysis 



approach. 



LAT PSR 


n 1 


^35 


/ 


/ 


T 


E 


B 






(10~ 8 cm" 2 s" 1 ) 


(Hz) 


(-10~ 12 Hzs" 1 ) 


(kyr) 


(lO^ergs- 1 ) 


(10 12 G) 


J0007+7303 


1509 


6.14(27) 


3.1658891845(5) 


3.6133(3) 


13.9 


45.2 


10.8 


J0357+32 


294 


0.64(10) 


2.251723430(1) 


0.0610(9)f 


585.0 


0.5 


2.3 


J0633+0632 


648 


1.60(17) 


3.3625440117(3) 


0.8992(2) 


59.3 


11.9 


4.9 


J1418-6058 


3160 


5.42(38) 


9.0440257591(8) 


13.8687(5) 


10.3 


495.2 


4.4 


J1459-60 


1089 


1.26(21) 


9.694596648(2) 


2.401(1) 


64.0 


91.9 


1.6 


J1732-31 


2843 


3.89(33) 


5.087952372(2) 


0.677(1) 


120.0 


13.6 


2.3 


J1741-2054 


889 


1.31(17) 


2.417211371(1) 


0.098(7) 


392.1 


0.9 


2.7 


J 1809-2332 


2606 


5.63(31) 


6.8125455291(4) 


1.5975(3) 


67.6 


43.0 


2.3 


J1813-1246 


1832 


2.79(24) 


20.802108713(5) 


7.615(4) 


43.3 


625.7 


0.9 


J 1826- 1256 


4102 


2.47(27) 


9.0726142968(4) 


9.9996(3) 


14.4 


358.2 


3.7 


J1836+5925 


2076 


8.36(31) 


5.7715516964(9) 


0.0508(6) 


1800.0 


1.2 


0.5 


J 1907+06 


2869 


3.74(29) 


9.378101746(2) 


7.682(1) 


19.4 


284.4 


3.1 


J 195 8+2846 


1355 


1.29(18) 


3.4436636006(8) 


2.6351(5) 


20.7 


35.8 


8.1 


J2021+4026 


4136 


10.60(40) 


3.769079109(1) 


0.7780(7) 


76.8 


11.6 


3.9 


J2032+4127 


2371 


3.07(26) 


6.9809351235(8) 


0.9560(4) 


115.8 


26.3 


1.7 


J2238+59 


811 


0.96(11) 


6.145017519(3) 


3.722(2) 


26.3 


90.3 


4.1 
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Frequency (Hz) 

Fig. 1. Frequency and frequency derivative distribution of the new pulsars. The dots represent the ^1800 pulsars 
in the ATNF catalog (20), the triangles represent the new pulsars reported in this paper, and the squares represent 
the 6 previously-known EGRET gamma-ray pulsars, including the Crab. The stars represent the new population 
of gamma-ray millisecond pulsars detected by Fermi-hPS and reported in the companion publication (10). The 
shaded region shows the parameter space covered by our blind search and includes ~ 86% of the ATNF pulsars. 
The three sets of lines illustrate the characteristic age r = —f/2f, inferred surface magnetic field strength B = 

3.2 x 10 19 ^/— /// 3 Gauss, and spin-down power given by E = A-K 2 Iff erg s~\ where I is the neutron star 
moment of inertia (/ = 10 45 g cm 2 ). This figure shows that the blind-search gamma-ray pulsars are drawn from 
the highest spin-down luminosity portion of the general population. 
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Fig. 2. Folded light curves, with a resolution of 32 phase bins per period, of the 16 pulsars discovered with Fermi- 
LAT, using 5 months of data with E > 300 MeV, selected from a region of radius 0.8° around the best position 
for the pulsar. The light curves are not background subtracted and can include a substantial contribution from the 
Galactic diffuse gamma-ray emission, particularly for the pulsars at low Galactic latitude. The x axis represents 
phase and the y axis counts per phase bin. Two rotations are shown and the phase of the first peak has been placed 
at ~0.3, for clarity. 
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